Frustrated spin chains in strong magnetic field: dilute two-component Bose gas regime 
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We study the ground state of frustrated spin-5 chains in a strong magnetic field in the immediate 
vicinity of saturation. In strongly frustrated chains, the magnon dispersion has two degenerate 
minima at inequivalent momenta ±Q, and just below the saturation field the system can be effec- 
tively represented as a dilute one-dimensional lattice gas of two species of bosons that correspond 
to magnons with momenta around ±Q. We present a theory of effective interactions in such a 
dilute magnon gas that allows us to make quantitative predictions for arbitrary values of the spin. 
With the help of this method, we are able to establish the magnetic phase diagram of frustrated 
chains close to saturation and study phase transitions between several nontrivial states, including a 
two-component Luttinger liquid, a vector chiral phase, and phases with bound magnons. We study 
those phase transitions immerically and find a good agreement with our analytical predictions. 

PACS numbers: 75.10.Jm, 75.30.Kz, 75.40.Mg, 67.85.Hj 



I. INTRODUCTION 

Frustrated spin systems have been an object of avid 
interest to researchers for several decades. In low dimen- 
sions, the interplay of frustration and quantum fluctua- 
tions proved to be favorable for generating states with 
unconventional spin order, such as chiral, nematic, or 
general multipolar phases, as well as for the existence 
of various disordered states. A strong external magnetic 
field competes with the exchange interaction, inducing 
a rich variety of phase transitions. An overview of this 
rapidly progressing field may be found in Refs. [l| and[^ 

A spin chain with competing nearest and next-nearest 
neighbor exchange interaction represents a paradigmatic 
model of frustrated spin systems in one dimension. De- 
spite many years of extensive studies, it continues to de- 
liver fresh surprises. The frustrated spin-i* chain model, 
which will be the main subject of our study, is described 
by the following Hamiltonian: 

n 

+ + "^^1+2) + J2^n^n+2^^ (1) 

where {S^,S^) are spin-S* operators acting at site n of a 
one-dimensional lattice, Ji, Jf and J2, J| are nearest- 
neighbor (NN) and next-nearest neighbor (NNN) ex- 
change interactions, and H is the external magnetic field. 
In this paper, we are interested in the frustrated case, so 
J2 is chosen to be positive and the sign of Ji is arbitrary. 
It is convenient to use the quantity 

P = J1/J2 (2) 

as the frustration parameter. The system may be alterna- 
tively viewed as two antiferromagnetic chains connected 
by ferro- or antiferromagnetic zigzag couplings Ji, Jf . 



For S* = i , the above model has been extensively stud- 
ied. Its isotropic (Jf = Ji, J| = J2) version has a rich 
magnetic phase diagram exhibiting states with compet- 
ing unconventional orders, both for a ferromagnetic'^"— 
and an antiferromagnetic^"^^ NN exchange interaction 
Ji. In particular, the vector chirality (VC), which is 
a quantum remnant of the classical helical spin order 
and is equivalent to the local spin current, competes, 
on the one hand, with multipolar orders characteriz- 
ing quasi-condensates consisting of multimagnon bound 
states, and, on the other hand, with the two-component 
Tomonaga-Luttinger liquid (TLL2)^^ where the two com- 
ponents correspond to magnons with momenta around 
two degenerate dispersion minima at inequivalent points 
in the Brillouin zone. The same set of phases has been 
found in the anisotropic version of the model as well.~ 

The model ([T]) with S — ^ and a ferromagnetic (FM) 
NN coupling has been discussed recently in connection 
with the description of several quasi one-dimensional 
magnetic materials such as LiCuV04 (Ref. [isl - fisl ). 
Rb2Cu2Mo30i2 (Refjli), Li2ZrCu04 (Ref. [IS), and an- 
hydrous CuCl2 (Ref. [21^. 

For higher values of S, not much is known about the 
magnetic phase diagram of frustrated chains. A theo- 
retical approach based on bosonizationS predicts that in 
the regime of two weakly coupled antiferromagnetic (AF) 
chains (|/3| ^ 1) a VC phase completely fills the magnetic 
phase diagram for all nonvanishing values of the magne- 
tization between zero and saturation. On the FM side 
(/3 < 0), for spin-S' > 1 chains right below the satura- 
tion field, a recent large- S" analysis in the framework of 
a dilute magnon gas approacbS^ predicts the existence of 
either a VC phase, for (icr{S) < /3 < 0, or a metamag- 
netic magnetization jump, for —4 < j3 < j3cr{S), with 
Pcr{S > 6) — —4. Numerical studies of magnetization 
curves^l have established the universal presence of a mag- 
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netization plateau at one third of the saturation value in 
antiferromagnetic chains with S = 1, |, and 2, while a 
metamagnetic jump at saturation has been found in FM 
chains^ for S = 1 and S — ^ (metamagnetic jumps have 
been observed earlier in 5 = ^ chains^ as well, particu- 
larly at saturation for (3 ^ —4). The numerical analy- 
sis of correlation functions^^'^'' has been limited to cuts 
at the fixed value of the frustration parameter /3 = 1, 
and has led to the conclusion that the field-induced VC 
phase is present in AF chains with higher S as well (val- 
ues up to S" = 2 have been studied). The TLL2 phase 
has been observed in 5 = ^ frustrated AF chains^iiir— , 
in anisotropic S — frustrated FM chain/'* as well as in 
S' = 1 chain with bilinear and biquadratic exchange.— 
Indications of its existence have also been found in the 
form of kinks in magnetization curves in higher-S" frus- 
trated AF chains. -'^ 

In this paper, we study the ground state of a strongly 
frustrated (|/3| < 4) spin-S" chain ([!]) in strong magnetic 
fields in the immediate vicinity of saturation. Just below 
the saturation field, the system can be represented as a 
dilute lattice gas of magnons, and since in the case of 
strong frustration the magnon dispersion has two degen- 
erate minima at inequivalent points ±Q in momentum 
space, this lattice gas is effectively a two-component one. 
We present a theory of effective interactions in such a 
dilute magnon gas that is especially adapted to the one- 
dimensional case, does not involve an 1/S expansion, and 
thus allows us to draw quantitative predictions for arbi- 
trary spin. We establish the magnetic phase diagram of 
isotropic and anisotropic frustrated chains close to satu- 
ration and study phase transitions between several non- 
trivial states, including TLL2 and VC phases, and phases 
with bound magnons. Particularly, we show that the 
TLL2 phase appears in FM and AF chains for all val- 
ues of the spin S. We complement our analytical results 
with density matrix renormalization group ^^i^^ (DMRG) 
simulations. Specifically, we compute the central charge 
as a function of /? (at fixed magnetization close to the 
saturation value) , which allows us to locate phase transi- 
tions between one- and two-component states. In order to 
identify those one-component states that have vector chi- 
ral order, we further calculate the chirality. Our DMRG 
results for the transition from the VC to the TLL2 phase 
are in a good agreement with our theoretical predictions. 

The structure of the paper is as follows: in Sect. [Tl] 
we describe mapping of the spin problem to the dilute 
two-component lattice Bose gas and present the theory 
allowing us to calculate effective interactions in such a 
gas for general spin S. In Sect. lIIIl the two-magnon scat- 
tering problem in a frustrated chain is considered, and 
its links to the effective theory of Sect. [ll] are established. 
Section|lV]discusses the specific predictions of the theory 
for frustrated spin chain models with different S and ei- 
ther isotropic or anisotropic exchange interactions, while 
Sect. |V] presents the results of numerical analysis and 
their comparison with analytical predictions. Finally, 
Sect. IVTl contains a brief summary. 



II. EFFECTIVE TWO-COMPONENT BOSE 
GAS MODEL 

We are interested in the interactions between magnons 
in the regime of a dilute magnon gas, i.e., for values of 
the field H just slightly below the saturation field Hg. 
Because of the diluteness of the gas, two-body interac- 
tions dominate. It is convenient to use the Dyson-Maleev 
representation^^ for the spin operators in ([Ij : 



s„ 



2Sbl 1 



2S 



(3) 



where 5„ are bosonic operators acting at site n. To en- 
force the constraint < 25", one can add the infinite 
interaction term to the Hamiltonian, which reads: 



n^n + uY^:{bib^) 



2S+1 



u 



(4) 



where : (...) : denotes normal ordering. Obviously, at 
the level of two-body interactions this term is important 
only ioT S — ^. 

Passing further to the momentum representation for 
bosonic operators, one can rewrite the model ([1]) in the 
following form: 

k kk' q 

(5) 

Here L is the total number of lattice sites, the chemical 
potential is /i = Hs — H , where 



25J2(l + /3V8) + 2S'(Jf + J|) 



is the saturation field^^ for the case of a strong NNN 
interaction | Ji | < 4 J2 , which is of our main interest here. 
The magnon energy 

Ek = 2S{Ji cos k + J2 cos 2k - Jl - J|) + Hs (6) 

is defined in a way that ensures that min£fe = 0. The lat- 
tice constant has been set to unity. It is easy to convince 
oneself that Sk does not depend on 2 ■ The dispersion 
Ek has two degenerate minima at fc = iQ, where 

Q = arccos [ - /3/4] . 

The two-body interaction Vq{k, fc') generally depends on 
the transferred momentum q as well as on the incoming 
momenta fc, fc': 

yffc fc'i-/'^'('^)~^['^(^) + '^(^')] m 



where 



J(fc) = 2 Ji cos fc + 2J2 cos 2fc, 
(fc) = 2Jl cos fc -I- 2 J| cos 2fc. 



(8) 
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— — 
k'-p 



The macroscopic couphngs F^^ are given by 



FIG. 1. The ladder approximation to the Bethe-Salpeter 
equation for the renormalized two-body interaction vertex 
Tq{k, k' ; E). Solid lines denote bare propagators. The ap- 
proximation becomes exact at /i 0. 



It is easy to see that in this notation Sk — S[j{k) — J{Q)\ . 

The model ([5]) defines a gas of bosons with a nontriv- 
ial double-minima dispersion, on a one-dimensional (Id) 
lattice. We are interested in the renormalized two-body 
interaction in this gas in the dilute limit, i.e., — > 0. 
In the — >■ limit, the self-energy vanishes,— so the 
full propagator coincides with the bare one. Then, the 
Bethe-Salpeter (BS) equation for the renormalized two- 
body interaction vertex Fg(/c, fc'; E) (where E is the total 
energy of the incoming particles) in the limit — > takes 
the following simple formi^ 



Tq{k,k';E)^V,ik,k') 



(9) 



1 ^ Vg-p{k +p,k'- p)Fp(fc, k'; E) 
T. 2^ 



E 



This expression is schematically shown in terms of Feyn- 
man diagrams in Fig. ([1]). For the physical "on-shell" 
vertex the energy is fixed at E ^ Sk + Sk' , but we keep 
Eq. ([9]) in a general form for reasons that will become 
clear shortly. 

At a low density of magnons (small chemical potential 
/i) the system is mainly populated by particles with mo- 
menta around the two dispersion minima at ±Q, which 
at low energies can be interpreted as two different bosonic 
flavors. One can obtain the low-energy effective theory in 
the form of the Gross-Pitaevsky-type energy functional 
for two-component field (for the outline of the derivation, 
see Appendix 1X1) : 



•Hgp 



dx 



E 

a=1.2 



o|2 



1, 



2m 



T^A 



(10) 



Here <i>^'^ are the macroscopic bosonic fields that describe 
magnons with momenta k lying within the intervals | k ± 
Ql < A around the dispersion minima, A is an infrared 
cutoff, and Ua — 1$"^ are the corresponding densities. 
Both bosonic species have the same effective mass 



k=Qj 



5J2(16-/?2 



(11) 
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r^(Q,g;0), 

r;}(-g,Q;0) + F^Q(-Q,Q;0) (12) 



The focus of our study here is to extend this 
For the sake of clarity, the Planck constant is set to unity. approach to the two-component gas, which, as we have 



where the vertex function T^{k,k' ; E) is the solution of 
the BS equation ([SJ with the infrared cutoff \p\ > A in- 
troduced into the sum over transferred momentum p. 

As noticed first by Batyev and Braginskii,^^ for a 
model with short-range interaction the solution of © 
can be expressed in terms of a finite number of Fourier 
harmonics in the transferred momentum. In our case, 
from the structure of Vq{k,k') it is easy to see that 
Tq{k, k'; E) generally contains harmonics proportional to 
1, cosg, sing, cos2q, and sin2q. The integral equation 
thus is reduced to a system of five linear equations that 
can be solved for general spin S*, thus no l/S* expansion 
is necessary. For the purpose of finding only the effective 
couplings Fii, F12, the number of linear equations in the 
general problem outlined above can be reduced from five 
to three, see Appendix [Bl 

The point Fn = F12 corresponds to the enhanced 
SU{2) symmetry at the level of the effective low-energy 
theory. For F12 < Fn, the ground state of the gas con- 
tains an equal density of the two particle species, and 
for F12 > Fii just one of the two species is present in 
the ground state. Recall that in our spin problem the 
total number of each species is not separately fixed, in 
contrast to a typical setup for atomic Bose mixtures. In 
a setup with fixed particle numbers the ground state at 
F12 < Fii is in the mixed phase, and F12 > Fn cor- 
responds to phase separation. In the spin language, in 
two or three dimensions the state with only one bosonic 
species corresponds to the helical magnetic order, while 
the state with an equal density of the two species cor- 
responds to the "fan" phase with coplanar spin order.— 
In one dimension, no long-range order is possible, and 
the state with only one species can be identified with the 
vector chiral (VC) phase, ^ while the two-species state 
corresponds to the two-component Tomonaga-Luttinger 
(TLL2) liquid. i3'25 

The approach that we have just described is well- 
known in the physics of magnetism in a strong field^i^— 
and is rather straightforward for three-dimensional sys- 
tems where in the limit A — )■ the couplings F^^ are 
well defined and acquire finite values. In the case of 
low dimensions, however, for A — >■ the integral in 
Eq. Q diverges at p = for (fc,fc') = {Q,Q), and for 
{k, k') = {—Q, Q) there are divergences both at p = 
and p = 2Q. Thus, the above scheme needs a modifica- 
tion in the low-dimensional case. The couplings F^j, can 
be interpreted as functions of the running cutoff A in the 
spirit of the renormalization group (RG) approach. The 
RG flow is then interrupted at a certain scale A = A* 
that depends on the particle density. This approach has 
proved to be quite successful for one-component Bose 
gas 
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(a) 




FIG. 2. The typical behavior of the running coupling con- 
stant r as a function of the infrared cutoff A, in a continuum 
model of a Boss gas with a bare contact interaction g (i.e., 
the two-body interaction potential is approximated by the 
delta-function U{x) = gS{x)): (a) for the case of the repulsive 
coupling g > 0, (b) for the case of attraction <? < 0. 



seen, is of special interest for frustrated systems. In this 
paper, we concentrate on the one-dimensional case. 

Solving the BS equation one can show that the 
expansion of F^f, as a series in A has the following form: 



1 


m 


1 


^ 11 


ttA 


.911 


1 


m 


1 


^ 12 


ttA 


912 



-0(A) 
-0(A) 



(13) 



Note that the expansion starts with the term propor- 
tional to 1/A that turns out to be the same for both 
couplings. Thus, for A ^ both F^^/A and F^j/^ Aow 
to the same universal value, which reflects the tendency 
of the RG flow to restore the SU{2) symmetry for the 
two-component Bose mixture.'*^ 

The parameters gn, gi2, which determine the second 
term in the expansions (jl3p . have a special meaning: un- 
der certain conditions they can be identified (see Appen- 
dices IAIC|) with the effective bare coupling constants of 
the continuum two-component Bose gas with a contact 
interaction 9iiS{x — x') between the same species (here 
5{x) is the Dirac delta- function) and an interaction of 
the form gi25{x — x') between different species. Such a 
continuum model is essentially the two-component Lieb- 
Liniger model. The behavior of running couplings F^j, in 
this problem is exactly given by the first two terms in 
Eq. ([T^."^^ In the semiclassical limit S — ^ oo, the param- 
eters gab are given by simple expressions^- 
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912 
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-2{1 + Ji + J^)^H,/S, 

J|) + 4(J| + 1), 



1 , 



(14) 



Typical behavior of F^^^ for the Lieb-Liniger model is 
shown in Fig. [2l It is clear that Tab changes its behavior 
on the characteristic scale Aq ~ gab^n. For the repulsive 
case {gab > 0) this is simply the scale at which the renor- 
malized coupling F^^ starts to deviate considerably from 



its bare value, while for the case of attraction {gab < 0) 
the renormalized coupling has a pole at A = Ag that 
technically limits the RG flow and physically indicates 
the existence of bound states with energy Eb ^ Aq/2to. 

It is important to realize the following: since we actu- 
ally deal with a lattice problem that has a natural ultravi- 
olet cutoff Auv ~ I7 the interpretation of gab as physical 
bare couplings makes sense only if the condition 



|gii|m<l, |gi2|m< 1 



(15) 



is satisfied (see Appendix [C] for more details). If the 
above condition of weak coupling is broken, parameters 
gab cannot be interpreted as physical bare couplings, and 
only the fully renormalized interactions 

i afc — i afc 

retain their meaning as effective low-energy coupling con- 
stants. For instance, if gab < but jgablm > 1, there is 
no attraction in the lattice problem and no bound states 
are present: directly looking at the behavior of F^j, for 
the lattice problem, one can convince oneself that there 
is no pole in A. In the strong coupling regime l^afcl'Ti > 1 
the sign of gab merely determines the curvature of the F^^ 
dependence near A = 0, but the renormalized coupling 
itself remains positive. The only physical meaning of gab 
in such a case is that they are connected to the asymp- 
totic phase shift of scattering states at small transferred 
momenta, see Sect. IIIll below. 

The scale A*, at which the RG fiow has to be inter- 
rupted, can be easily fixed from the well-known exact 
result for the single-component dilute Bose gas:44 in the 
limit /Lt — >■ the particle density n must behave as 



n = {2m^i/^T^f'^ 



(16) 



Indeed, in the phase with only one bosonic species, which 
corresponds to F12 > Fn, the density of particles ob- 
tained from the semiclassical description of Eq. (fTO]) is 
n = /u/Fii, so in order to fuUfil Eq. one has to 

assume th at at /i — > the effective coupling behaves as 
Fii — > TTyJ ^i/{2m) and thus 



A, — \/ ^m/2. 



(17) 



In the phase with equal densities of different species, 
which corresponds to F12 < Fn, the total density at 
the semiclassical level of Eq. (|10p is n = 2/i/(Fii -I- F12). 
It follows from Eq. (|13p that for /i — )■ the behavior 
of the density is described by the same expression (fT6|) 
as for the one-component gas. Thus, at the transition 
from the vector chiral phase to the two-component TL 
liquid, in the regime of a low magnon density there is no 
magnetization kink (cusp), in agreement with numerical 
results j^^ii^ The reason for that is the effective fermion- 
ization of the magnon gas and formation of a single Fermi 
sea at low densities, see the discussion below in Sect. IIIII 
From Eqs. (IT5|) one can see that the transition point 
between the TLL2 and VC phases, that is determined by 
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the condition Fn = T12, can be detected by the crossing 
of gi2 and gu (i.e., a change of sign of gi2 — gn that may 
happen if one of the couphngs goes through a pole does 
not correspond to a phase transition). We will consider 
specific examples below in Sect. IIVI 

There is another way of calculating the renormal- 
ized interaction in low dimensions, which avoids intro- 
ducing the infrared cutoff. In Refs. 45 and ^[6, it has 
been shown that a good approximation to the correct 
"many-body" expression for the effective four-point ver- 
tex function (i.e., the two-particle scattering amplitude 
in the presence of a finite particle density) at low en- 
ergy r*^^(fc, fc'; E^O) can be obtained from the "bare" 
two-particle scattering amplitude Tq{k, k'; E) taken at a 
finite negative energy E = —E^ that is proportional to 
the chemical potential fi: 

rf^(fc,fc';i;->o)wr,(fc,fc'; -£;,), E^^Cfi. (i8) 

Here the coefficient C turns out to depend on the dimen- 
sion d of the problem: for d = 2 one obtains C = 1, and 
C = 7r^/8 for d = \. The two-particle scattering ampli- 
tude Tq{k, k'\ E) is precisely the solution of Eq. ^ since 
/i = means that no other particles are present except 
for the two that scatter off each other. 

Following this scheme, we can define the energy- 
dependent coupling constants as 

T^^{E)=Vo{Q.Q\E), 

yi2{E) = ro(-Q, g; E) + r2Q(-Q, Q; E) (19) 

(note that the expressions on the rhs of Eq. ([T5)) are 
now calculated without any infrared cutoff). Solving 
the BS equation in the symmetrized form (jB2[) . (|B3I) for 
^ ab{—E.t), one can show that the expansion of r^j,^(— £'*) 
in E^ has the following structure: 
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Vi2{-E.) \4EJ gi2 ' ^ ' ' ^^"^ 

The similarity with the expansion in the infrared cutoff, 
i.e., Eq. (IT5|) . is obvious. Replacing i?* with its Id value 
Eif = 7r^^/8, one can readily convince oneself that 



ah 



Tab{-E,) + 0{ii- 



3/2N 



(21) 



i.e., both regularization methods yield equivalent results 
in the dilute limit — >■ 0. However, for practical pur- 
poses, the off-shell method is more convenient since it 
allows an easier calculation of all integrals involved. 



III. TWO-MAGNON PROBLEM IN 
FRUSTRATED SPIN CHAINS 

A general two-magnon state with the total quasi- 
momentum K can be written in the form 



I* 



K, 



= J2 ae^^("+^ 



rt_ rt- 



\F) 



n,r>0 



{E\S^^^Sn Sn Sj^^j,\F) 



where \F) is the fully polarized state. The Schrodinger 
equation for the two-magnon problem leads to a coupled 
system of equations for the amplitudes Cr, presented in 
Appendix IdJ that is a generalization of the one for the 
case of an SU (2) symmetric exchange. Both scattering 
states and bound states of two magnons can be obtained 
from the solution of this system. The bound states can be 
obtained along the lines derived for an SU (2) symmetric 
exchange Here we discuss the solutions for the case of 
continuum states that describe scattering of two magnons 
with momenta fci, ^2- 

Solving this scattering problem, we have to distinguish 
two cases. In the first case the problem is non-degenerate, 
meaning that for a given total momentum K ~ ki + ^2 
and energy E there is only one scattering state. This 
will be realized when the momenta of the two magnons 
participating in the scattering are in the vicinity of the 
same dispersion minimum, e.g., fci = Q + k, and k2 = 
Q — k. In this non-degenerate case the system (jD3p can 
be solved by using the following ansatz: 



Cr 



; (fcr -I- (5ii) -I- ve 



r > 1. 



(22) 



This ansatz satisfies the recurrence relation in Eqs. (jD3|) . 
for (?■ > 3) if Ko is a root of the following equation: 

K 

Ji cos — (cos /c— coshKo)+'^2 COS K{co& 2fc— cosh 2kq) — 0. 

(23) 

Furthermore, to be a physically acceptable solution, the 
correct root should satisfy Re[Ko] > 0. The three un- 
known constants Co, w, and 5ii are determined from the 
first three equations of (|D3p . 

To extract the scattering length an from the phase 
shift one can use the same relation as found in the 
one-dimensional continuum model of a one-component 
Bose gas interacting via a short-range potential, namely 



<^ii(fc) 



— flii/c at fc — > 0, or, in a different form. 



an — lim \ cotiSii) I k\ . 



(24) 



If the short-range potential between two particles of mass 
TO with coordinates x, x' is approximated by the contact 
potential of the form U{x — x') — gii6(x — x'), then the 
scattering length an is connected to the bare coupling 
gii as follows >^ 



gu = -'2/iaiim). 



(25) 



In the second case, there is more than one scatter- 
ing state for a given energy and total momentum. This 
degeneracy stems from the parity symmetry and the 
double- minimum structure of the single-magnon disper- 
sion. This case will be realized if the momenta of two 
scattering magnons are in the vicinity of two different 
dispersion minima. For instance, the state with the to- 
tal momentum if = that corresponds to fei ^ Q + p, 
k2 = —Q — p (half relative momentum k = Q + p) is de- 
generate with another scattering state with ki = Q — p, 



FIG. 3. (Color online). Single magnon and two-magnon 
spectra for the isotropic spin-1 chain with /? — —2.5 at the 
field H > Hs, as a function of the total momentum K. The 
solid line shows the single-magnon dispersion with minima at 
±Q. One can clearly see degeneracies of the scattering states 
in the vicinity of the lower bound of the two-particle contin- 
uum for K — as well as in some high-energy regions. A 
peculiar feature of the spectrum is the existence of a stable 
two-magnon bound state inside the continuum of scattering 
states (this happens for any spin S > ^ and /3 < 0). The 
procedure of finding two-magnon bound states in isotropic 
ferromagnetic S > | frustrated chains was discussed in detail 
in Ref. [22, and here we include them for the sake of complete- 
ness. 

k2 = —Q + p (half relative momentum k = Q — p)- 
This degeneracy is visible at the lower bound of the two- 
magnon scattering continuum around if ~ 0, depicted 
in Fig. [3] Generally, p ^ p because the shape of the 
one-magnon dispersion curve is not exactly symmetric 
around the minima at ±Q, but p — > p at p — 0. The 
relative momenta k, k of two degenerate states with the 
same total momentum K are connected by the equation 

K ~ 
Ji cos — (cos k — cos fc) -I- J2 cos A'(cos 2k — cos 2k) — 0. 

(26) 

For K = the solution of the above equation takes the 
form 

k = arccos(2 cos Q — cos k), (27) 

from which, in particular, it follows that for p = (k — 
Q) ^ one ha.s p = {Q - k) = p + 0{p^). 

The Schrodinger equation for the scattering states 
Eq. (jD3|) is in this case solved by the following ansatz: 

Cr = cos (fcr -I- a{k)) -f v{k) sin(fcr), r > 1. (28) 

From the form of the scattering state Eq. (l28l) it follows 
that if one prepares the incoming state of two magnons at 
t — > — 00 with relative momentum fc, the outgoing state 
at t ^ -\-oo will be a superposition of states with relative 
momenta k and k. 

To determine the three unknown quantities Co,a, and 
V we again use the first three equations of Eq. (jD3[) . In 
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FIG. 4. (Color online). Bare coupling constants for S = 
f/2: solid lines correspond to the calculation based on the 
regularized Bethe-Salpeter equation, while symbols show the 
results extracted from the solution of two-magnon scattering 
problem. In the region jS > 0.614, the intraspecies coupling 
(711 enters the super- Tonks regime. The interspecies coupling 
g\2 enters the super- Tonks regime in the region P > —2. 

Appendix[E]we show how to extract the interspecies scat- 
tering phase shift 

r cos a 

012 — — arccos _ (29) 

a/ cos'^ a + [v -\- sin a)^ 

from the scattering state Eq. (PH)) . The scattering length 
is then obtained in analogy to the case of the intraspecies 
scattering, 

ai2 = lim [cot((5i2)/_p] (30) 

and the interspecies coupling is given by a formula com- 
pletely similar to Eq. (PS)) : 

5i2 = -2/(TOai2). (31) 

It is worth noting that scattering states vanish when 
p — > 0, i.e., when the magnons' momenta tend to the 
dispersion minima fci Q, k2 ^ ~Q. This is a con- 
sequence of fermionization at low densities and indicates 
that different "species" in fact form a single Fermi sea 
at low energies. The existence of a single Fermi sea is 
the reason for the absence of a magnetization cusp at the 
transition between the vector chiral and TLL2 phases4^ 
It is worth noting that a similar phenomenon is known in 
the one-dimensional Hubbard model: in the limit of low 
electron density the electrons with different spin form a 
common Fermi sea 4^ 

The resulting coupling constants for S = ^ chain are 
shown in Fig. S] One can see that on the antiferromag- 
netic side (Ji > 0) the scattering length an is typi- 
cally small (of the order of the lattice constant) and goes 
through zero (i.e., gn has a simple pole). Similarly, on 
the ferromagnetic side there is a pole in 512. In the previ- 
ous section, we have seen that if the condition Eq. (|15l) is 
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not satisfied, a mapping to the effective continuum model 
of the Lieb-Liniger type (with contact interactions) is 
problematic. In relation to that, we would like to discuss 
the scattering length for an unfrustrated ( J2 = J| =0) 
spin-i XX Z chain, which shows a similar behavior. Set- 
ting Ji = 1, one can write the analytical expression for 
this scattering length as 

axxz = Ji7(l + Jl) ■ (32) 

Close to the ferromagnetic point Jf — —1, it shows res- 
onant behavior and this allows one to develop the effec- 
tive Lieb-Liniger model as the continuum limit in this 
regioui^ At Jf = 0, however, axxz goes through zero 
and this can be interpreted as follows: the continuum 
limit of the XX Z spin-i chain for Jf > does not fit 
into the description in terms of the repulsive Lieb-Liniger 
model, rather, it is a super Tonks gas (sTG)^ with the 
Luttinger parameter K < 1. Including now the second 
neighbor interaction J| (but not the hopping J2), one 
obtains that at Jf — ^ cxj and J| —?> 00, the scattering 
length goes to 2. Thus, it seems that in the regime of 
strong coupling |(7|m ^ 1 the scattering length has a dif- 
ferent physical meaning, namely it can be interpreted as 
the excluded volume of a particle. What is the correct 
continuum limit in this case, however, remains a compli- 
cated issue. 



IV. ANALYTICAL PREDICTIONS FOR 
FRUSTRATED SPIN-S CHAINS 

The formalism developed in Sect. |lT] is rather general 
and applies to any lattice Bose gas with degenerate in- 
equivalent dispersion minima. Let us look in more detail 
at the predictions of our theory for the frustrated spin 
chain model defined by the Haniiltonian ([1]). 

Figures [1][5] show the values of the "bare couplings" 
gii and gi2 for frustrated spin-S" chains with isotropic 
interactions (i.e., Jf = Ji, Ji — J2) with S = ^, 1, and 
I as functions of the frustration parameter /3 = Ji/J2- 
One can see that even for the lowest spin values S — ^ 
and 1, there is an excellent agreement between the ap- 
proach based on the regularized Bethe-Salpeter equation 
and the results extracted from the direct solution of the 
two-magnon scattering problem. 

A. Ferromagnetic frustrated isotropic chains 

(Ji < 0, A = 1) 

On the ferromagnetic side {/3 < 0) the intraspecies cou- 
pling gii remains small enough for its interpretation as 
the bare coupling of the effective continuum Lieb-Liniger 
model to be justified (strictly speaking, this is true for 
5 > 1 where |(7ii|to <C 1, while for S" = i the coupling 
values are higher, |(7ii|m < 0.5). For S" = i, this cou- 
pling remains negative in the entire range of —4 < /3 < 0, 
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FIG. 5. (Color online). Bare coupling constants for isotropic 
frustrated spin chains with (a) S' = 1 and (b) S = \: solid 
lines correspond to the calculation based on the regularized 
Bethe-Salpeter equation, while symbols show the results ex- 
tracted from the solution of two-magnon scattering problem. 



which indicates the presence of bound states and is con- 
sistent with known numerical results i^ii^ 

For 5 = 1, the intraspecies coupling 51 1 becomes neg- 
ative in the region — 4 < /3 < /?/ ~ —2.11 (see Fig. [5^), 
which signals the appearance of bound magnon states. 
As shown previously in Ref. [l^ by means of an l/S" ex- 
pansion and numerical analysis, this also leads to meta- 
magnetism (a finite jump in the magnetization at the 
saturation field). The finite value of this jump can be 
understood in our approach by noticing the fact that the 
full renormalized effective coupling F^f depends on the 
particle density n — 1 — M, where M is the magnetization 
in units of saturation, and for [3 < f3f the coupling F^-j* 
becomes positive at M below some critical value which 
determines the magnitude of the jump in M that is nec- 
essary to make the system stable again. The values of 
M that are spanned by the jump correspond to states 
with magnons collapsed into a single "drop" . Such states 
have negative magnetic susceptibility and are avoided in 
the "grand canonical" setup when one fixes the external 
magnetic field H . 

The numerical results of Ref. Il^ also suggest that the 
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FIG. 6. (Color online). The critical value Sct of the spin 
that marks the disappearance of the metamagnetic phase, as 
a function of the order n of expansion in The exact 

value is 5cr = 5, while the leading-term calculation presented 
in Ref. [S^ gives Scr ~ 6. 

collapsed state has a finite vector cliirality. This is easy 
to understand by observing that the interspecies inter- 
action remains repulsive in the entire ferromagnetic re- 
gion — 4 < /? < 0. Thus, it is energetically favorable to 
form a "drop" involving only one magnon species. As 
a side note, in this case gn can be directly interpreted 
as the bare coupling of the effective continuum theory, 
even though the interspecies coupling 1712 remains large 
I512IW > 1. This can be done because the system is 
effectively in a one-component state. 

Frustrated ferromagnetic chains with spins up to S* = | 
behave quahtatively in the same way as described above 
for 5 = 1. The interspecies "bare coupling" 512 goes 
through a simple pole only for S* = i (for 5 = 1 it goes 
through the pole of order 2 at /3 = 0), so for S* > 1 it 
becomes positive everywhere (see Figs. H1 |5I7|) . The in- 
traspecies coupling gn remains small and vanishes at 
/3 = /?/. Upon increasing 5', the transition point /?/ 
quickly moves very close to ^ = —4 and disappears for 
S > 5. This is clear from the behavior of 1711 at /3 — > —4: 

911 = + 4)^ + ((/? + 4)^/2) . (33) 

Thus, Scr = 5 is the critical value of spin where the meta- 
magnetic behavior vanishes in isotropic chains (for S > 5, 
it exists only in the presence of an easy-axis anisotropy) . 
We would like to note that Ref. [H reported a slightly 
different value of 5'cr = 6; this discrepancy is due to the 
fact that Ref. [13 used just the leading term in the large- 
S expansion, while our present approach is exact to all 
orders in 1/S. Figure [6] illustrates the behavior of ^cr as 
a function of the order of the l/S expansion. 

B. Antiferromagnetic frustrated isotropic chains 

(Ji > 0, A = 1) 

On the antiferromagnetic side < /? < 4, the promi- 
nent feature of isotropic chains is that for realistic (low) 
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P=J/J, 

FIG. 7. (Color online). Typical behavior of the "bare cou- 
plings" gii and (J12 for isotropic spin chains with small and 
large spin S as functions of the frustration parameter /3, com- 
pared to the semiclassical, S — 00 result (shown with thin 
solid line for gn and thin dashed line for (?i2): (a) for small 
S, there is a single crossing of and gi2 curves, while (b) for 
large S > 12, those curves cross three times (see the inset). 

values of the spin S the condition Eq. (|15l) is not satisfied, 
as can be seen from Figs.|3][S] Moreover, gu changes its 
sign through a pole in this interval of /3, for any value 
of S. As discussed above, in this regime the system can- 
not be mapped to an effective continuum Lieb-Liniger- 
type model; however, a mapping to the effective classical 
Gross-Pitaevskii model (fTOl) is still possible, and the cor- 
responding effective couplings Fn and F12 remain posi- 
tive in the entire antiferromagnetic region. For large 5*, 
the mass m cx S~^ becomes small (except near the two 
Lifshitz points (3 ~ ±4) which ensures that in the limit 
S ^ 00 the condition Eq. ([15]) is again satisfied every- 
where except for small regions around the pole and the 
antiferromagnetic Lifshitz point. 

Figure [7] shows the typical behavior of the "bare cou- 
plings" gii and gi2 for isotropic spin chains with small 
and large spin S, compared to the semiclassical S* — ?> 00 
result. In the 5 — > 00 limit, gi2 > gu everywhere, and 
512 = 511 only at two points /3 = —4 (the ferromag- 
netic Lifshitz point) and /3 = 2 (the Majumdar- Ghosh 
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FIG. 8. (Color online). Phases ol isotropic frustrated antifer- 
romagnetic (Ji > 0) spin-S* chains just below the saturation 
field, as predicted by our effective two-component Bose gas 
theory. The region /3 > 4 corresponds to a one-component 
TL liquid phase (TLLI). If one varies /3 at fixed S, there 
is a single transition between the two-component TL liquid 
(TLL2) and the vector chiral (VC) phase for S < 12, while 
for S > 12 one encounters three consecutive TLL2-VC transi- 
tions. Symbols correspond to the transition points extracted 
from our numerical analysis of chirality correlators, see Sec- 
tion |V] for details. 
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1/2 


1 


3/2 


2 


12 


/^TLL2-VC 


1.682 


1.688 


1.694 


1.702 


1.80, 2.89, 3.13 



TABLE I. Critical values of the exchange coupling /3 — J1/J2 
for the transition from the VC to the TLL2 phase in isotropic 
frustrated antiferromagnetic chains of spin 5* at the saturation 
field, from the dilute Bose gas theory of the present work (cf. 
Figure |8j . 
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FIG. 9. (Color online). Bare coupling constants for S — ^ 
anisotropic chain with / Ji = JI / J2 = A: (a) J1/J2 = —0.3 
and (b) J1/J2 = -3.0. 



point) 43 For small S, there is only one crossing of gu 
and gi2 at some /? = /3i < 2, indicating a transition be- 
tween the VC phase that exists in the interval < /? < /?i 
and the TLL2 phase that occupies the region /3i < /3 < 4. 
With increasing S, the transition point /3i moves toward 
the Majumdar-Ghosh point /3 = 2, see Fig.|8]and TableU) 
For very large S > 12, two more crossings of the gu and 
gi2 curves appear (see the inset in Fig.[7l)), and another 
piece of the VC phase emerges in the region 2 < /3 < 4, 
filling it completely as S goes to infinity (see Fig. |8]). 

C. Anisotropic chains (A 7^ 1) 

Motivated by the availability of numerical results^'' 
for anisotropic frustrated ferromagnetic chains with S = 
1/2, we would like to check our theoretical predictions for 
anisotropic systems as well. Figure |9] shows the behav- 
ior of the "bare couplings" for S* = ^ anisotropic chains 
with two different values of /3 < 0, as functions of the 
anisotropy Jf/Ji = JI/J2 = A. For /3 = -0.3, the 
intraspecies coupling gu remains small and vanishes at 
A ~ 0.56, while gi2 is negative and large (I512I™ ^ 1). 
This indicates a transition from the vector chiral phase at 



A < 0.56 to a phase with bound states at A > 0.56. This 
agrees favorably with the numerical results of Ref. [13 
that reports a transition between the chiral and nematic 
phases at A ~ 0.6. 

For P = —3, there is a crossing of gn and 512 at A ~ 
0.5 (see Fig. I^b) that indicates a transition between the 
TLL2 phase at A < 0.5 and the VC phase at A > 0.5, 
nicely fitting to the numerical data of Ref. [13. Both gu 
and gi2 go through poles and become positive (which 
does not change the phase since in the course of this 
change one always has ri2 > Fn), then they become 
small and positive around A = 1, and gu crosses zero at 
A ~ 0.947, indicating a transition from the VC phase to 
a phase with bound states (312 changes its sign later, at 
about A = 1.05). Again, this is in good agreement with 
the numerical study Ref. [l3 that observes a transition 
around A = 0.95. 

Figure [TO] shows the predicted phase diagram of 
anisotropic S — ^ and 5* = 1 chains immediately be- 
low the saturation field in the {13, A) plane. One can see 
that with increasing S the TLL2 phase shrinks and the 
phase with bound states is pushed toward higher A (for 
5 > 5 it lies entirely in the easy-axis region A > 1). The 
vector chiral phase dominates the phase diagram. One 
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FIG. 10. (Color online). Phase diagrams of anisotropic 
(Ji/Ji — J2/J2 ~ A) frustrated chains with (a.) S — ^ and 
(b) S = 1 right below the saturation field, derived from our 
dilute Bose gas approach. The boundary between the two- 
component Tomonaga-Luttinger liquid (TLL2) and the vec- 
tor chiral (VC) phases has been extracted from the condition 
gii = gi2. At the boundary of the region labeled "Bound 
states" the intraspecies coupling gn changes its sign from pos- 
itive to negative (it turns out that gn always becomes nega- 
tive before 1712 does) . The "bound states" region for 5 = | in 
fact contains multiple phases corresponding to bound states of 
a different number of magnons,— which is beyond the scope 
of the present theory that can only detect the presence of 
two-magnon bound states. For S = 1, the "bound states" 
region is metamagnetic,— i.e., it exhibits a finite magnetiza- 
tion jump immediately below the saturation field. For S = ^, 
circles correspond to the numerical data of Ref. [T3 . and the 
diamond is the location of the VC-TLL2 transition at A = 1 
as found numerically in Refs. [ll| and [T^. For S = 1, sym- 
bols denote our DMRG results (see Section |V] for details): 
diamonds and circles denote the phase transition points ex- 
tracted from the analysis of the chirality correlators and the 
magnetization curves, respectively. 



can conclude that the present theory reliably detects the 
transition points in anisotropic ferromagnetic chains even 
at such low values of the spin a.s S = ^. It is worth 
noting that the transition point between the TLL2 and 
VC phases in the (3 = —3, S = ^ chain comes out cor- 
rectly despite the fact that both |(7i2|m- and |5ii|to are 
large at the transition. On the antiferromagnetic side, 



the only numerical data available in the literature is for 
the isotropic S — ^ chain. ^^'^^ In this case, for a reason 
that is not yet clear, the agreement between the theory 
and numerics is considerably worse. 



V. NUMERICAL ANALYSIS 

To verify our theoretical predictions, we have per- 
formed extensive density matrix renormalization group 
(DMRG) simulations. For a detailed description of the 
DMRG technique?^ we refer the reader to the review 
Ref.iai We focus on two quantities: First, we calculate 
the central charge from the scaling behavior of the von 
Neumann entropy, which allows us to distinguish one- 
from two-component phases. Second, we compute the 
vector chirality order parameter from the analysis of chi- 
ral correlation functions, to identify those one-component 
phases that have long-range vector chiral order. The goal 
is to provide a numerical check of the phase diagrams 
for isotropic and anisotropic frustrated chains. Figs. [8] 
and 1101 respectively. In order to identify the region with 
bound states in anisotropic chains, we calculate magneti- 
zation curves. We have studied systems of different sizes 
L ranging from L — 50 to L — 400, both with periodic 
and open boundary conditions, typically keeping 600 to 
1500 states in most calculations. 



A. Central charge 

Our theory suggests that close to the saturation field, 
in antiferromagnetic frustrated chains with spin values of 
1/2 < 5* < 12, one encounters two phase transitions when 
increasing /3 from zero: (i) from the VC phase into the 
TLL2 phase and (ii) from the TLL2 phase into the TLLl 
phase (compare Fig. [S]). The TLL2 phase has central 
charge c — 2, while the VC and the TLLl phases have 
c = 1, suggesting that a calculation of c can be helpful 
in determining the phase boundaries. 

To accomplish that, we exploit that the von Neumann 
entropy S^n can be obtained as a by-product of any 
DMRG calculation. Furthermore, from the finite-size 
scaling of S^n, one can extract the central charge. The 
von-Neumann entropy is defined as 



SvN,l{1) = -tr(/9; Inpi) , 



(34) 



where pi is the reduced density matrix for a subsystem 
of length I embedded in a chain of a finite length L. In 
a gapless state that is conformally invariant, the I and L 
dependence of the von Neumann entropy is given 



SvN,l{1) = g In 



'L . 




— sm 







'5, 



(35) 



which is valid for systems with periodic boundary con- 
ditions (PBC). Here g is a non-universal constant that 
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FIG. 11. (Color online). The von Neumann entropy Si, jv,l(0 
vs block length Hor 5 = 2 at M = 7/8 for /3 = 2 (circles) and 
P — 1 (squares). The solid lines are fits to Eq. (|35p . resulting 
in c « 1.01 and c — 2.10, respectively (the last 10 data points 
were excluded from the fits). Symbols are DMRG data for 
m — 1500 states. 



depends on the magnetization M. As DMRG di- 
rectly accesses the eigenvalues of these reduced density 
matrices r^iSl it is straightforward to measure SyN,L/2{l) 
with this method. In the case of open boundary condi- 
tions (OBC), the von Neumann entropy is given by one 
half of the formula Eq. ([55)1 . where besides g, there are 
also contribution that oscillates as a function of I. 

Typical DMRG results for the von Neumann entropy 
of chains with L — 64 spins and PBC are presented in 
Fig. [m for 5 = 2, together with a fit of Eq. to the 
numerical data (lines). These fits result in c « 1 for 
^ = 1 and c sa 2.1 for ^ = 2 at M = 7/8. As it turns out, 
there are strong finite-size effects in the estimation of c 
from fits to Eq. (|35p . in particular at high fields close to 
saturation which are our primary interest. Nishimoto^ 
has recently suggested to calculate c from the following 
formula: 



ln[cos(7r/L)] 



(36) 



which turns out to converge much faster with the system 
size. However, the numerical data needs to be sufficiently 
accurate for Eq. (|36| to yield reliable results, and we 
found that one has to push the number of states in a 
typical DMRG calculation to about m ~ 1500 to achieve 
good accuracy. For instance, for the parameter of Fig.lTTl 
Eq. dSe]) resuhs in c = 1.005 for ;3 = 1. We used Eq. §6^ 
to extract c whenever the DMRG calculation proved to 
be sufficiently accurate, and resorted to fits to Eq. ([55)) 
(or its OBC counterpart) otherwise. 
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FIG. 12. (Color online). Typical vector chirality correlation 
functions for an isotropic S — 1 chain of L = 200 spins. 



B. Chirality order parameter 

To study the behavior of the vector chirality, we have 
calculated correlators of the form 

(KiKj), Ki = [Si X Si+lY ■ (37) 

The chiral order parameter k has been extracted as 

= lim {KiKj)av, (38) 

where the index "av" means that correlators were aver- 
aged over pairs of sites (i, j) with the same |« — j|, while 
care was taken to stay "in the bulk", i.e., i and j were 
kept sufficiently far away from the ends of the chain. Fig- 
ure[T2]shows typical vector chirality correlation functions 
for a S" = 1 isotropic frustrated chain of length L — 200. 

C. Numerical results: Isotropic chains 

Our results for the central charge of isotropic frustrated 
chains are shown in Fig. [131 We have collected data at 
fixed magnetization Af for S" = 1, |, 2. From the figure 
we clearly identify two phase transitions: (i) at a critical, 
S'-dependent /? > 3 from TLLl (c = 1) phase into the 
TLL2 phase (c = 2) and (n) at /3 > 1.6, from the TLL2 
phase into another c = 1 phase. Based on our theory, we 
expect this latter phase to have vector chiral order (which 
is verified by calculating the chiral order parameter, see 
below) . 

The position of the phase transition from TLL2 into 
VC depends only weakly on S, consistent with our the- 
ory (compare the values given in Tab. [Tj). According to 
our theory, the phase transition from TLLl into TLL2 
happens at approximately /3 = 4, the Lifshitz point, in 
the limit of M ^ 1. Our numerical data for M < 1, 
however, shows that this transition sets in at a slightly 
smaller /3 < 3.5 than what one would expect from Fig. [5) 
This suggests that the phase boundary depends on M. 
Indeed, for instance, in the case of 5 = 1, this phase 
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FIG. 13. (Color online). Central charge at M = 7/8 for 
5 = 1, 3/2, and 2, as extracted from Eq. (|36l) [The data 
points in the c ~ 2 region, i.e., 1.5 < /3 < 2.5 were extracted 
from direct fits of Eq. ^ to the DMRG results for 5'„jv,l(0]- 
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FIG. 14. (Color online). Square of the vector chirality order 
parameter for frustrated isotropic antiferromagnetic S — I 
chain at fixed value of the magnetization M — 95/100, as a 
function of the frustration parameter /3 = J1/J2. 



boundary shifts to smaller values of /3 as M decreases 
which can be deduced from the analysis of magnetization 
curves presented in Ref. H^, where this particular phase 
transition shows up as a kink in M{H). We have not 
systematically studied the M dependence as the DMRG 
data for SyN,L tends to converge rather slowly at large 
M . Our numerical results for c are in very good overall 
agreement with the theory. 

We next verify that the c — \ phase at /3 < 1.6 is the 
VC phase, by directly calculating the chiral order param- 
eter as described above; the results are shown in Figs.fTH 
and[T51 One can clearly see from Fig. [15] that, in agree- 
ment with our analytical predictions, the transition point 
P from the vector chiral phase to the two-component 
Tomonaga-Luttinger liquid increases slowly with the in- 
creasing spin value S. Our numerical results for the tran- 
sition boundaries in isotropic frustrated chains are shown 
as symbols in the phase diagram of Fig. [51 

Figure [TJj illustrates that with the increasing system 
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FIG. 15. (Color online). Square of the vector chirality order 
parameter for frustrated isotropic antiferromagnetic spin- 
s' chains, at fixed value of the magnetization M = 9/10, as 
a function of the frustration parameter f) = J1/J2: (a) for 
chains consisting of L = 100 spins, at different values of S; 
(b) for S = 12 chains of different length L. 

size L, the VC-TLL2 transition becomes more steep. 
This is consistent with the first-order character of this 
transition as predicted by the theory (recall that the VC- 
TLL2 transition corresponds to phase separation in the 
language of the effective two-component Bose gas the- 
ory). 

For high magnetization sectors it is possible to study 
large spins as well and results for S" = 12 are included in 
Fig. [121 One should take those S = \2 results with some 
caution: our data for small systems (consisting of 80, 100, 
and 120 sites) still show strong nontrivial size dependence 
(see Fig. [T5b) that precludes infinite-size extrapolation, 
and we could not study larger systems because of the 
rapidly growing dimension of the Hilbert space. Nev- 
ertheless, the resulting picture for S" = 12 qualitatively 
agrees with our analytical prediction, namely, one can see 
the recurring TLL2-VC phase transitions in accordance 
with the phase diagram of Fig. [51 

D. Numerical results: Anisotropic 5 = 1 chain 

We have studied numerically the phase diagram of 
the anisotropic (JI/J2 = Jl/Ji — A) frustrated spin- 
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FIG. 16. Square of the vector chirality order parameter 
for frustrated antiferromagnetic S = 1 chain for A = 0.7, at 
fixed value of the magnetization M = 9/10, as a function of 
the frustration parameter /? = Ji/J2- 



1 chain, to verify the analytically predicted phase dia- 
gram of Fig. [TUb . Figure [1^] shows the typical depen- 
dence of the chiral order parameter k on the frustra- 
tion parameter /3 for anisotropic 5=1 chain, at fixed 
anisotropy A and magnetization Af , for three different 
system sizes. One can notice that finite-size effects are 
not always monotonous, which complicates the extrapo- 
lation to infinite system size. 

In Fig.[l7l we present our results for the dependence of 
the central charge on j3, for the same parameters as used 
for the chirality in Fig. [161 These data were obtained 
from systems with open boundary conditions for which 
we could obtain reliable data. The additional oscillatory 
contributions to S^n are evident from the inset of Fig. 
[TBI The behavior of the central charge is consistent with 
the results for the chirality: Whenever > 0, c = 1. 
At both /3 > 2 and /3 < —3.5, where k drops to zero, 
the central charge goes up to c = 2, indicating phase 
transition from VC to TLL2 phase. 

As for the regime of almost decoupled chains |/3| <C 
1, one should keep in mind that this is a particularly 
difficult parameter region for DMRG calculations, and 
the apparent deviation from the analytical predictions 
(according to which c = 2 is only realized for one point 
/3 = 0) is an artefact of the rapidly deteriorating accuracy 
at /3 ^ 0. 

In order to study the phase boundary between the VC 
and "bound states" (metamagnetic) phase, we have also 
calculated magnetization curves, and looked at the value 
AS^ of the steps by which the magnetization changes. 
While inside the VC phase (where magnons interact re- 
pulsively) the magnetization grows in steps of AS*^ = 1 
when increasing the magnetic field, after crossing the 
metamagnetic region boundary the height of steps in the 
magnetization curves quickly changes to higher values 
AS*^ — 2, 3, etc. Typical results for M{H) are shown 
in Fig. [iSl In this example [S — 1, (3 = —1-2), we see 
metamagnetic behavior (a jump in M{H)) directly un- 
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|3=J/J, 

FIG. 17. Main panel: Central charge, obtained with open 
boundary conditions, for a frustrated antiferromagnetic S = 1 
chain for A = 0.7, at magnetization M = 9/10, as a function 
of the frustration parameter /3 = J1/J2. Inset: von Neumann 
entropy SvN,l(1) versus i at /3 = 1,2.6. 
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FIG. 18. Typical magnetization curves for S = 1 chain 
with P = —1.2, at different values of the anisotropy A. For 
A = 1.1, the curve exhibits a metamagnetic jump at the 
saturation field (similar jump for A = 1, for —4 < j3 < —2.11 
has been already discussed in Ref.— ). 



der saturation magnetization for A = 1.1, but not for 
A = 1.02 or A = 0.98. The results for the onset of 
the metamagnetic behavior as a function of A and j3 are 
included in Fig. [TOT b) (circles). 

Our numerical results for the TLL2-VC phase bound- 
ary are shown by symbols in Fig.fTOb. The overall agree- 
ment between the theory and numerics is very good. 



VI. SUMMARY 

We have developed the effective theory for low-density 
one-dimensional lattice Bose gas with two degenerate dis- 
persion minima. This study has been motivated by the 
problem of the ground state phase transitions in frus- 
trated spin-S* chains subject to strong magnetic fields just 
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slightly lower than the saturation field, where such a sys- 
tem can be viewed as a dilute lattice gas of magnons 
in the fully polarized background. At low density of 
magnons they can be treated as bosonic particles, and 
in a wide region of parameters the magnon dispersion 
has two degenerate minima at inequivalent points zbQ in 
the momentum space. 

The low-energy effective theory has the form of the 
Gross-Pitaevskii functional for a two-component field. 
We have developed the scheme that allows the calcu- 
lation of renormalized effective interactions Fn, ri2. 
Applying this formalism to frustrated spin- 5* chains, 
we establish the magnetic phase diagram of isotropic 
and anisotropic, ferro- and antiferromagnetic frustrated 
chains close to saturation and study phase transitions 
between several nontrivial states, including the two- 
component Tomonaga-Luttinger liquid (TLL2), vector 
chiral (VC) phase, and phases with bound magnons. Our 
scheme does not involve an l/S expansion, and thus al- 
lows us to draw quantitative predictions for arbitrary 
spin. 

We complement our analytical predictions by a vari- 
ety of numerical results obtained by means of the density 
matrix renormalization group technique. Particularly, we 
compute the central charge as a function of the frustra- 
tion parameter close to saturation which allows us to lo- 
cate phase transitions between one- and two-component 
states. Further, we directly compute chirality correlators 
in order to identify those one-component states that have 
vector chiral order. In addition, we study the magneti- 
zation curves to detect the presence of magnon bound 
states. These numerical results are in a good quantita- 
tive agreement with our theoretical predictions for chains 
of spin S > 1. At the same time, in the case of frustrated 
antiferromagnetic S = ^ chains, the quantitative discrep- 
ancy between our theoretical prediction and the numer- 
ical results by Okunishi^^ and Hikihara et oIf^ for the 
location of the TLL2/VC phase boundary is unexpect- 
edly large. Nevertheless, even for S ~ \ our theory still 
correctly captures the topology of the phase diagram. 
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Appendix A: Low-energy effective theory and 
connection to the two-component Lieb-Liniger model 

To obtain the low-energy effective theory for the 
Hamiltonian ([SJ, we define the bosonic field V' = 



S/c ^fcC**^^ with the Euclidean action given by = 
J dT{J dx'4>*dr'4' + "^[V"]}: and divide the field V into 
"slow" parts and the "fast" part 0: 

ip^Y^ 'J'"e""°'3^ (Al) 

a=l,2 

where A is the running cutoff. Treating the "slow" part 
as the background field, one can integrate over the fast 
modes / 2?[0]e^-^''' = e~-^* and obtain the effective ac- 
tion that is essentially the generating functional of 
vertex functions^ and can be written as follows: 

_|_ id,".* ^P^* tM):°'P 

^ 2^k+q,n+u^k'-q,n'-u:'~^q,k,k' ,n+n'^k,n^k' ,n' 
+ (higher-order terms in $). 

In the above equation, the summation over repeated in- 
dices is implied, and the A-dependent coefficients [/'■"•' 
can be expressed through the vertex functions F*^"^ of 
the initial theory determined by the action A,p'. 

='5„^r(2)(a,Q + fc,c.), (A3) 
C^S^'Tcf = rW(a„g + k, a^Q + k'; w), 
f^Sfc"? = r^KQ + k^'^.Q + k';oj) 

+ AZ~<^c)Q+k'~k^qi'^c.Q + k, apQ + fc';a;), 

with the restriction that sums in the internal lines of 
the corresponding Fcynman diagrams include only "fast" 
modes with momenta p outside the regions |p ± Qj < A. 

The two-point function — F*^^^ (fc, cj) is equal to the 
Green function G(fc,w). In the limit /i — > 0, i.e., 
in the absence of particles, the self-energy vanishes— 
and the full propagator G(fc,w) is given by its bare 
value G'^^-' (fc, oj) = (iw — Sk)'^- The four-point vertex 

Fq''^(fc, fc',a;) = Tq{k,k';E = iuj) is determined by the 
Bethe-Salpeter equation, which in the limit /i — )• takes 
the simple form of Eq. with the infrared cutoff A 
around the momenta ±Q introduced into the integral. 

We neglect higher-order terms in the expansion (jA2p . 
which is justified at low densities, replace e±Q+fc that 
enters C/'-^-' by its leading small-fc term fc^/(2m), and 
neglect the momentum and energy dependence of U^'^\ 
which is appropriate in the limit of small A. Then the 
effective low-energy action (|A2p takes the form of the 
Gross-Pitaevskii functional for a two-component field: 

A^^ j dr j dx^^°'*{dr- ^iW 

^|9^$"p-Ki$"*$'^*F;^^$"$''}, (A4) 

where the renormalized running couplings are given by 
F^o — C^o^o'^Q. As shown in the main text by solving the 



15 



BS equation, in the one-dimensional case the dependence r^(fc, fc'; E) = r_^(fc', k; E), one readily obtains the fol- 
lowing equation for Fg^^-*: 



of those couplings on the cutoff has the form 
= irA/m + Co^pK^ + ..., A ^ 0. 



(A5) 



Now, consider the system of two bosonic species Lpi^2 
interacting via contact potential, described by the action 



(A6) 



A^^ j dr j dx^ip*^{dT ~ y)ipa 

+ I^gal3na{x)ni3{x)^, 



where = |<PqP, a — 1,2 are the densities of the 
two components. This is nothing but the two-component 
Lieb-Liniger model. The low-energy effective theory for 
this model can be again cast in the form (jA4p .^^ with the 
renormalized couplings given by 



afi 



Ja0 



ttA 1 



ttA 



A->0. 
(A7) 

Comparing (jA7p and (|A5|) . one can establish a corre- 
spondence between the original lattice problem ([5]) and 
the simplified continuum model (jA6l) of the Lieb-Liniger 
type, that is characterized by bare couplings ga/s- 



n"'w=>-„(o.o)-ii: '^""5""'^' > 

^11 (<?, P)^l [V,-p{Q +P,Q-p) + iP^ -P)] , (B2) 



(12). 



and a similar equation for F, 



Ff)(£;) = 2Fi2(g,g)-yE 



VMq,p)r^p'\E) 
2e„ + E 



Vi2{q,p) = 2 [Vg-p{p, -p) + V,+p{~p, +p)] . (B3) 
Explicitly, the symmetrized kernels for our model are 



Viiiq,p) 

Vi2{q,p) 

where 



V^^((7,p) + [/iCosp + C/2Cos2p, S>1 
V''{q,p) + U, U^oo S'=i 

'V^q,p)-J{p), S>1 
V'{q,p) + U, U^oo 5 = i 



(B4) 



Appendix B: Symmetrized equations for Fn, F12 

To simplify the equation (j9|), we first symmetrize the 
kernel. Let us introduce the following functions that are 
even in the transferred momentum q: 

F(")(£;) = F,(Q,g;£;), (Bl) 
T(^^\E) = Fq+,(-Q, Q; E) + Fq_,(-Q, Q; E), 

then Tii{E) ^ F^"^ and Ti2{E) = fJLq- From (H]), us- 
ing the identities eu — S-k, Vq{k,k') = V-q{k' ,k), and 

I 



Ui = -2 Ji cos Q, U2 = -2J2 cos 2Q, (B5) 
V''{q,p) = 2Jf cos(g) cos(p) + 2J| cos(29) cos(2p) 



Solutions of Eqs. (IB2p . (IB3|) can be now sought in the 
form containing only even Fourier harmonics: 



r(") = A1+A2 cos q + As cos 2q, 
r(i2) = B1 + B2 cos q + Bs cos 2q. 



This ansatz transforms the integral equations (jB2l) . (jB3l) into 3x3 linear systems: 



25'—! 
2S{Ui+U2)-E 

~12 



'22 ^ 2jf 
~23 

^11 + 25(;7i+;72)-B 12 



'13 



'22 2,/^^ 



~13 








~23 


A2 






^33 + 2J| / 








~13 






^ 2/cs 


~23 


B2 




2 cos Q 


^33 + 2J| / 


\BsJ 




y2 cos 2Q 



where cs — I — 



E 



2S{Ui + U2) 

4e 



V ^13 

for 5 > 1 and C5 = 1 for S = ^, and the matrix coefficients are given by 



fifj 



£Q+k + £Q-k — E 



fifj 



L^2ek-E' 



/l = 1, /2=COSfc, /3 = COs2fc 



(B6) 



(B7) 



(B8) 



(B9) 



r 



The solutions of those equations can be easily obtained in a closed though somewhat cumbersome form. 
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Appendix C: The continuum limit for a 
one-component lattice Bose gas: bare and 
renormalized interaction 



Consider Bose gas on a lattice, with the single-particle 
kinetic energy Sk that behaves quadratically at small mo- 
menta, Sk — k'^ /2m at fc — > 0. Let us assume there 
is an on-site interaction of strength qq. In the limit 
of vanishing density, one can write down the Bethe- 
Salpeter equation for the renormalized two-body inter- 
action rq(fc = 0, fc' = 0; i? = 0) = r in the same way 
as we have done in Sect. HIl but now it is considerably 
simplified since the bare interaction vertex Vg(fc, k) = go 
is independent of momenta, and the renormalized inter- 
action r is independent of the transferred momentum q: 



r-5o(i-r/), / = 



1 



dp 
2^ 



(CI) 



Here Ayv = tt is the natural ultraviolet lattice cutoff, and 
A is the running infrared cutoff. The flow eventually gets 
stopped at A = A* given in Eq. pT| . At low gas density 
n the regulator A* = 7rn/2 <C 1, so we are interested in 
the behavior of r(A) at A -> 0. 

For the continuum Bose gas model with contact inter- 
action of the form U{x — x') = gc5{x — x') (the Licb- 
Liniger model), the equation for the renormalized inter- 
action Tc is very similar. 



J. r T\ T m f dp m 



(C2) 



The running coupling in the continuum model is thus 
given by 



9c 



(C3) 



Let us come back to the lattice expression (IC1[) . Since 
Sk — > k"^ /2m at fc — > 0, it is easy to see that the integral 
/ in (ICll) behaves as 



m 
ttA 



A-> 0, 



(C4) 



where /lat is some parameter of the order of unity that 
incorporates all effects of the lattice (ultraviolet cutoff 
as well as the deviation of the dispersion Sp from the 
quadratic law). The important point is that this con- 
stant /lat, depending on the lattice details, can be nega- 
tive: for instance, if one takes the free particle dispersion 
Sp ~ p^/(2m) and simply introduces the lattice cutoff 
Auv, then /lat — — m/(7rAuv)- Curiously, this constant 
/lat does not appear if the particle dispersion is a pure 
cosine (i.e., only nearest-neighbor hopping is taken into 



account): it is easy to show that for = ^(1 — cosp) 
the integral / in (|Cip is 



/ 



(1-^ + ...) at A^O, (C5) 



27rtan(A/2)) ttA' 
i.e., in this case /lat — ^ 0. 

The resulting behavior of the running coupling in the 
lattice problem can be represented in the form identical 
to that of the continuum model: 



F = 



30 



1 + 



r5o 



50 = 



30 



1 + 5o/i 



lat 



(C6) 



Comparing equations (jC3l) and (|C6p . one can see that go 
plays the role of the effective bare coupling gc if the orig- 
inal lattice problem is mapped to the continuum Lieb- 
Liniger model. The value of % will, e.g., determine the 
phase shifts in two-body scattering. As long as the cou- 
pling constant g^ is small, such a mapping presents no 
problem. However, it is easy to see that if /lat is nega- 
tive and fl'ol-'^iatl becomes close or exceeds 1, then even the 
sign of 5o can differ from that of the lattice bare coupling 
50- The parameter g^ can be directly interpreted as the 
effective continuum bare coupling only if the condition 



\9o\m < 1 



(C7) 



is satisfied, otherwise one has to look at the dressed cou- 
pling F(A,). 



Appendix D: Schrodinger's equation for the 
two-magnon scattering problem 
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We set the magnetic field at the saturation value 
H Hs = 2S{Jf + Ji) + S{Jf + 8J|)/4J2 and con- 
sider a general two-magnon state that can be written in 
the following forni (here we use a normalization different 
from that of Ref.lH): 



I* 



2M 



n.m>n 



I 



{F\SmSt Sn Sm\F) 



(Dl) 



where \F) is the fully polarized state. The total quasi- 
momentum K is conserved, therefore one can sepa- 
rate the center of mass motion by setting Cn.m = 
exp{i^^^^^Y^K)Cm-n- Solving the Schrodinger equation 
'H|^2Af) = E\'i'2M) for scattering states in the thermo- 
dynamic limit, one can identify the energy E — e{ki) + 
£(fc2) + Ef, where Ep = NiS^{Jf + Ji) - HsS) is the 
energy of the fully polarized state, as 

E = no-2SJ{Q) + EF, 
K 

n^/AS = Ji cos(A:) cos + J2 cos(2fc) cos(/i:), (D2) 

where ki, k2 are the magnon momenta, k — ^{ki — k2) 
is one half of the relative momentum, and K — ki + k2 
is the total momentum. The Schrodinger equation leads 
to the following linear system for the amplitudes C^: 
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QqCo = 2y^S{2S - l){JiCi cos y + J2C2 cos if), 



K 



K 



(17o - Jl)Ci = 2JiyJS{2S -l)Co cos y + 2S'JiC2 cos — + 2S'J2(Ci + C3) cos 



(17o - J|)C2 = 2J2V'5'(2S'- l)Co cos a: + 2S'J2C4 cos AT + 2SJi{Ci + C3) cos 

VloCr = 2SJi cos y (a+i + Cr-i) + 2S'J2 cos K {Cr+2 + Cr-2) > 3). 



(D3) 



r 



Appendix E: Scattering phase shift in the 
interspecies scattering problem 

To extract the phase shift for the interspecies scatter- 
ing problem from the solution (|28p. we put k = Q+p and 
consider the limit of small deviation from the minimum 
of dispersion p ^ 0, then k = Q — p + 0{p'^) and the 
scattering state Eg . ([25)) can be rewritten as follows, 



Cr— \J cos^ a + (w + sin a)^ cos(Qr) cos(pr + 6\2) 



— \J cos^ a + (w — sin a)^ sin (Qr) sin(pr — 7) (El) 



where 



(5l2 = 

7 = 



cos a. 



arccos 



arccos ■ 



•y/ cos^ a + (u + sin a)2 ' 
cos a 

•\/ cos^ a + (u — sin q:)^ 



(E2) 



From the form (jETj) one can easily guess that for Ji = 
(Q = 7r/2) the scattering phase shift is 5\2- This is 
indeed so, because to obtain the phase shift we have to 
put r = 2n and the second term in Eq. (|E1|) vanishes 
(for odd values r = 2n + 1 there is no scattering involved 



since at Ji = the system corresponds to two decoupled 
chains) . 

For (5 = (i.e., at the ferromagnetic Lifshitz point 
/3 — —4) and Q — tt (the antiferromagnetic Lifshitz point 
/3 — A) it is also easy to recognize S12 as the scattering 
phase shift. 

Furthermore, for general values of /3 one can show that 
limp^o ct = ±7r/2 and limp_^o v = ±1, so that limp_j.o(i' — 
sin a) = limp^ocosa = 0. This corresponds to the fact 
that scattering states vanish when p — 0, i.e., when the 
magnon momenta tend to the dispersion minima /ci — >■ Q, 
/c2 —7- — Q, and reflects the ferminionization of excita- 
tions forming a single Fermi sea, as mentioned in the 
main text. It is easy to see that the second amplitude in 
Eq. (jEip tends to zero, limp_>o \/ cos^ a + {v — sin a)^ — 
0, while the first amplitude tends to a constant, 
limp^o \/ cos^ a + {v + sin a)^ = 2. Thus, the scattering 
phase shift for interspecies scattering at p — > is given 
by di2 for any /?. 

One can also show that the two following expressions 
are equal to each other. 



lim 



cot {5i2) 



In other words, if a 
S12 -7r/2 + ai2p. 



= lim 



cot (a) 
2p 



— 7r/2 + 2ai2P at p 



(E3) 
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